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Abstract 



', The projection postulate has been used to predict a slow-down of the time evolu- 

Q | tion of the state of a system under rapidly repeated measurements, and ultimately a 

freezing of the state. To test this so-called quantum Zeno effect an experiment was 
performed by Itano et al. (Phys. Rev. A 41, 2295 (1990)) in which an atomic-level 
measurement was realized by means of a short laser pulse. The relevance of the 
results has given rise to controversies in the literature. In particular the projection 
■ postulate and its applicability in this experiment have been cast into doubt. In this 

paper we show analytically that for a wide range of parameters such a short laser 
' pulse acts as an effective level measurement to which the usual projection postu- 

late applies with high accuracy. The corrections to the ideal reductions and their 
accumulation over n pulses are calculated. Our conclusion is that the projection 
postulate is an excellent pragmatic tool for a quick and simple understanding of the 
slow-down of time evolution in experiments of this type. However, corrections have 
to be included, and an actual freezing does not seem possible because of the finite 
£h ' duration of measurements. 

El PACS numbers 03.65.Bz; 42.50.-p; 32.80.-t 

> : 

^ ■ 1. Introduction 

The so-called quantum Zeno effect (QZE) [|TJ is a theoretical prediction for the behavior 
of a system under rapidly repeated measurements at times At apart. It is based on usual 
quantum theory and on the concept of instantaneous measurements together with ensuing 
state reductions according to the projection postulate of von Neumann and Liiders @. 
The predictions are: 

1. Impediment and slow-down of the time development of the system due to repeated 
measurements. 

2. Freezing of the state for At — > 0, i.e. in the limit of continuous measurements. 

The underlying reason for this can be traced to the fact that for short enough times 
transition probabilities grow quadratically with time. If, in a given time interval T, one 
performs n measurements at times At = T/n apart, then the probability to find an 
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orthonormal state is at most proportional to n(At) 2 = T 2 /n, which goes to zero for 
n — > oo or At — > 0. 

Properties 1 and 2 may be taken as a definition of the QZE, as for instance in Refs. 
@ 0- & • A slightly different definition is used in Ref . M where essentially only property 
2 is used. 

An experiment to test the QZE for atomic systems has been performed by Itano et al. 
||, following a suggestion by Cook [^]. In the experiment a large number iV of ions was 
stored in a Penning trap (see Fig. 1 for the relevant level structure, a V configuration; 
level 2 is (meta-) stable). The time development is given by a so-called ir pulse M of 
length T n , tuned to the 1-2 transition frequency. A very short pulse of a probe laser 
couples level 1 with an auxiliary third level, and this is regarded as a measurement to 
which the projection postulate is applied as follows. It is assumed that an atom is either 
in level 1 or level 2, depending on whether or not it has emitted photons during the 
probe pulse ||- At the end of the n pulse the final level populations for up to 64 probe 
pulses ("measurements") during a 7r pulse (see Fig. 2) are determined and found to be in 
agreement with the quantum Zeno predictions fTOf . 

The interpretation of this experiment, as to whether it does or does not bear on the 
QZE, has been controversial in the literature. Some have hailed the results as a dramatic 
verification of the QZE, others argue that they are unrelated to it [§], ||, 11], 0, [13, 1||, 



15], [TBI, O. In particular the use of the projection postulate and reduction of the wave- 



function have been criticized, and the very existence of the QZE has been cast into doubt. 
It is no exaggeration to say that the QZE has aroused tremendous interest in the literature 
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It was pointed out in Refs. [§, ||, [TJ], [15[ [16|, |TjJ that these results could be understood 
without any recourse to the projection postulate or to the QZE. One can simply incorpo- 
rate the short probe pulses in the dynamics by an appropriate term in the Hamiltonian 
or in the optical Bloch equations. A numerical solution of these Bloch equations should 
then yield the experimental result, and indeed they do to good agreement. The projection 
postulate does not seem to be needed. Is it also incorrect to use it here? 

This is one of the main two questions we are going to address in this paper. First we 
give a justification of the projection postulate as a useful and approximate technical tool 
in experimental situations of the general setup considered in Ref. || for a wide range of 
parameters. Our motivation differs somewhat from that of Refs. @, [5|, 14, [15[ [If], [TTf] , 



although there is no contradiction, in particular not to Refs. 0, [T6j. In Ref. the 
ensemble density matrix is calculated by adiabatic Bloch equation techniques for a single 
probe pulse and found to be nearly diagonal. But this contains no information about the 
outcome for individual atoms and subensembles and does not prove that atoms with or 
without photon emissions are in |1) and |2), respectively. The same technique is used in 
Ref. |IJ| to study the parameter domain of "good" measurements with a resulting diagonal 
density matrix We, however, are concerned first with selective measurements, namely 
with the states of atoms at the end of a probe pulse with or without photon emissions. 
These states are explicitly determined for the first time and found to be close, but not 
identical, to 1 1) and |2) [^2j. Thus a probe pulse can indeed be regarded as accomplishing 
a highly accurate - but not perfect - reduction to and measurement of levels 1 and 2 
in experiments of the type of Ref. H for a wide range of parameters. As a second 
question we discuss the cumulative effect of the deviations from ideal measurements - i.e. 
measurements which can be described by the projection postulate - over n probe pulses 
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for the density matrix and determine the resulting level populations at the end of the tt 
pulse. We exhibit an explicit n dependence of the corrections to the treatment by the 
projection postulate. This is also a new result. 

For a probe pulse to be an effective measurement some obvious requirements have to 
be fulfilled. First of all, its duration, r p , should be very short compared to the duration 
T w of the it pulse, 

r p «T, . (1) 

Furthermore, the probe pulse cannot be too short or too weak, because it should produce 
fluorescence photons from level 1 with high certainty. In addition, at the end of a probe 
pulse any population in the auxiliary level 3 should decay to level 1 extremely rapidly. 
This means that A^ 1 , the inverse Einstein coefficient of level 3, must be tiny compared 
to the time between two probe pulses, 

A- z 1 « ^ - r p . (2) 

The first two conditions lead to a restriction on two parameters, namely [23] 

^3^2 , Q 2 / \ 

c p = ~nT~ ^ 1 and e R = FT ^ 1 • ( 3 ) 

"3 ^3 

In Sections 4 and 5 we will use the conditions 

e p = -j-^ < l a nd e d ee -f < 1 . (4) 

Then the condition 1 is automatically fulfilled since e R = ^e 1 / 2 . This parameter 



regime is compatible with the one investigated in Ref. |16[]. In the experiment of Itano et 
al. § e R is about 6.5 • 1(T 6 , e p about 4.1 • 10" 4 , and e d about 2.5 • KT 4 |2J. 

In our paper we mainly use the quantum jump approach (quantum trajectories) 



27], ^3] , and because of its inherent simplicity the analysis can be carried out analytically 
here for a wide range of experimental parameters. This approach is equivalent to the 
the Monte-Carlo wave function approach f29 |. The ensemble of all trajectories satisfies 



the usual Bloch equations and is often used to obtain numerical solutions of the latter 
through numerical simulations. The quantum jump approach deals with pure states for 
single atoms instead of density matrices for ensembles |[25|| , which in the present case means 
a simplification from nine components to three, and it describes the time development of 
single atoms between photon emissions by a simple reduced Hamiltonian and by a "jump" 
to a reset state (in this paper the ground state) at a photon detection [BO . It allows an 
intuitive understanding of the processes involved. 

The simple derivation of the quantum jump approach outlined in Section 2 uses the 
projection postulate and reductions as a technical tool for photon detections (not for 
atomic measurements!), and it might seem that one uses this postulate for an investiga- 



tion of itself. However, as explained elsewhere |26[ the use of reductions in the derivation 
of quantum trajectories is not necessary and can be avoided. Moreover, the photon detec- 
tions happen on a much shorter time scale than all times considered in the experiment, 
and for such detections the projection postulate has been a reliable tool in the past. 

In Section 3 we consider the simple case of a probe pulse with the n pulse temporarily 
switched off. In this example it is very easy to see how the probe pulse acts on a state 
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aii|l) + ct 2 1 2) of a single atom and that it effectively leaves it either in |1) or in |2), 
with probability |«i| 2 and |a2| 2 , respectively, depending on whether the atom has emitted 
photons or not. The small deviations from an ideal measurement with perfect state 
reduction are determined. 

In Section 4 we consider the case of a single probe pulse with the 7r pulse switched on. 
Complications arise now since the 7r pulse causes a small additional transition between 1 
and 2. In Section 5 the case of n probe pulses during a single n pulse and the build-up of 
the corrections are considered. It turns out that one gets an excellent approximation for 
the level populations if in the results for the idealized case one takes the finite duration 
of the probe pulse into account by neglecting the action of the n pulse during this time 
(see Fig. 3). This simply means that T^/n is replaced by AT = T n /n — r p . At the end 
of Section 5 we compare our analytical results with those of a numerical solution of the 
optical Bloch equations for the parameters of the experiment. The agreement is amazing. 

In the last section we discuss our result and their significance for the role of the 
projection postulate in the so-called quantum Zeno effect. Our main conclusion is that, 
although the projection postulate is not necessary here, it is a useful tool for a fairly 
accurate description of the measurements involved. It is useful since it gives a quick 
intuitive understanding of the physical situation - without having to solve unwieldy Bloch 
equation. Insofar the projection postulate should not be dismissed out of hand. However, 
it is only approximately valid, because in practice realistic measurements are never ideal 
nor instantaneous. Therefore the first part of the Zeno effect - impediment and slow-down 
of the time development - can legitimately be understood by the projection postulate, 
but the second part - the freezing of the state - is in our opinion an over-idealization. 



2. The quantum jump approach in quantum optics. Quantum trajectories 



In this section we briefly summarize the quantum jump approach used in the subsequent 
sections. The reader familiar with it can proceed directly to Eqs. ([16]) - fll9|). The idea 
is to describe the radiating atom between photon detections by a reduced (or effective) 
time evolution operator giving the time development under the condition that no photon 



has been detected ||26|| . After a photon detection one has to reset the atom to the reset 
state ("jump"), with ensuing reduced time development, and so on. For a driven system 
with many emissions one then obtains a stochastic path, called a quantum trajectory |28| . 
The general reset state has been determined in Ref. f2~7|| . For a V system as considered 
in this paper the reset state after an emission is the ground state. The reduced time 
development together with the reset state provide a complete stochastic description of 



the time development of the atom [27]. Starting with this description one can then derive 
the Bloch equations describing an ensemble of radiating atoms p6| , [27]]. In fact both 
approaches, quantum jumps and Bloch equations, are possible and equivalent ways to 
describe the time evolution of an ensemble of fluorescing atoms, but the former is also 
easy to apply to the emission behavior of a single atom. 

We now indicate how to determine the reduced time development operator for the 
V system. To be sure that no photon has been detected in a time interval one may 
imagine measurements on the radiation field in a rapid succession at times At apart, at 
t\ < ■ ■ ■ < t m = t say. If in all these measurements no photons are found the state of the 
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atom at time t is, by the von Neumann - Liiders projection postulate ||31|| , 

\<P(t)) = K h }(o ph \u(t m ,t m ^)\o ph }(o ph \... 

...K^iOphWi^to^OpM) (5) 

= \O ph )U ied (t,t )\iP). (6) 

where \0 p h) is the vacuum state, |0 P ^) the initial state, U the complete time development 
operator, and where the second equality serves as a definition of the reduced time evolution 
operator U Ted (t,t ), which acts on atomic states. The time difference between successive 
measurements has to be chosen short enough to be able to say that at most one photon 
has been detected in this interval. On the other hand this time difference has to be longer 
than the inverse transition frequencies [PEfl . Under these assumptions one can calculate 
Eq. (||) by means of perturbation theory. With U xe &{t, t ) one has a simple expression for 
the probability Po(t) that no photon is detected in the interval [0,t] if the atom is in the 
state \ip) at t — 0, 

Po(t) = |||0(t))|| 2 

= \\u Ted (t,t )m\ 2 ■ (7) 

We now apply this to the V system depicted in Fig. 1. In this system the upper levels 
2 and 3 couple to a common ground level 1, with Einstein coefficients A 2 and A s . Later 
we will consider level 2 to be stable and put A 2 = 0. We assume here that U32 = ^3 — 0^2 
is in the optical range, i.e. not too small. For simplicity we consider zero detunings of 
the driving fields, whose (real) Rabi frequencies are denoted by fl 2 and fl 3 , respectively. 
The Hamiltonian in rotating wave approximation is given by |32j 

H = H° A + H° F + j2 7i{~a|i)<l|e-^ + h.c.} + £ £ % {g ikX a kX \i) (1| + h.c.} 

i=2 1 i=2 kA 

where 

1 /2 

g ikX = i eD a ■ e kA (u; kA / (2e hV)) , (8) 

with the transition dipole moment Y>n = (z(X|l) and e k A the polarization vector. V is 
the quantization volume, later taken to infinity. Going over to an interaction picture with 
respect to 

H = H° A + H° F (9) 

one has 

Hi{t) =EM^K)(l|+h.c.} + ^^^{^ kA a kA K)(l|e i ^-^^ i + h.c.} . (10) 

By Uj we denote the corresponding time development operator. With At = ti — £j_i 
in the range given above we can calculate in second order perturbation theory the time 
evolution under the condition that no photon has been detected. From 

= \ A - % - t ^(Opftlffj^lOp^-i t dt' f dt"{% h \H I {t')H I {t")\% h ) 
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one obtains for the first-order contribution 

3 1 



-iY^{Mi)(l\+h.c.}At 

2 



and for the second order, omitting terms proportional to (At) 

rtf 

>ti— i Jti-i 



i,i=2 kA 



= "EE rate**-"'* [ t '~ U ~dTg^ x g; kX e^-^\z}(j\. (11) 

i,i=2 kA 

In the last equation we have substituted r = t' — t" '. Since At is much larger than the 
inverse optical frequencies cjj -1 one can extend the inner integral to infinity, leading to 
7r5(u)j — oukx) plus a principle value. Alternatively one can argue that the correlation 
function 

%(r)=E^* kA e^-^ )T (12) 

kA 

is negligible for r ^> uj~ . The sum over k then yields generalized decay constant and 
level shifts A^-. The level shifts are small |33|] and will be neglected in the following. With 

e 2 D,-i Di, , . . 

= -sAnrS (13) 



one obtains for the second-order contribution 

3 r At 



E 1 './ ')'•/ / dre~^-^ T . (14) 

i,j=2 J ° 

Now, if u 3 — u 2 is in the optical range, as supposed here, then the last integral vanishes 
for i j and equals At otherwise. We note that 

L i% 2 1 

where A\ is the Einstein coefficient of the i— th level. Collecting all terms we thus obtain 

3 1 3 1 

<O ph |LM*i,*i-i)|(W = ^-^{^^(il+WAt-EoAI^X^IAt. (15) 

i=2 Z i=2 Z 

This can be written as exp{— i H^ ed At/h} where, in matrix notation, the reduced Hamil- 
tonian H^ cd and the atomic operator M are defined through 



, o n 2 n 3 

HUn = -\£l 2 -iA 2 j =-%M. (16) 
" ; <h -iA 3 



Later on, we will take A 2 = 0. For arbitrary time intervals we thus have, in the interaction 
picture, 

C// ed (t, 0) = e ~ lH ^ t/h = e~ Mt . (17) 
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The no-photon probability is then, for t = and initial state ip, or more generally a 
density matrix p(0), 

P (t;V) = l|e- Mt |^)H 2 (18) 

P (*;p(0)) = tr{e- M V(0)e- Mti } . 

The probability that the first photon is emitted in (t, t+dt) equals Po(t; ip) —Po(t+dt; ip) = 
Wi(t; ip)dt where 

wi{t;i>) = -j t P (t;ip) (19) 

is the probability density for the first photon |54" |. For small upper level separation 
nonzero off-diagonal terms may appear which lead to interesting coherence effects 



35| , S7[ |38j . For general n-level systems the reduced Hamiltonian is given in Ref . |27| . 

The reduced time development is not unitary. The reason is that it does not describe 
the time evolution of the whole ensemble but that of the subensemble with no photons. 
The size of this subensemble is decreasing in time since an atom for which a photon has 
been detected leaves the sub-ensemble, and this is reflected by the decrease of the norm 
squared in Eq. (|7]). The above probability density determines the (random) time for the 
first photon. After that the atom is reset to the ground state, |1), for a V system. The 
next emission time is then determined by wi(t; 1), and so on. In this way one obtains a 
quantum trajectory. 

From this description of single systems one can recover the usual Bloch equations of 
the complete ensemble as follows f27|. The density matrix p(t) of the ensemble is a sum of 
two terms, p > and p°, corresponding to a subensemble of atoms with or without photon 
emissions until time t, respectively. From Eq. ([17]) one has 

p°(t;p(0)) = e- M V(0)e- Mt ' . (20) 

If /(t; p(0))dr denotes the (unconditioned) probability to find a photon between r and 
r + dr, then the sub-subensemble of atoms with their last emission before t in this interval 
is described by 

l(r;p(0))drp°(t-r;\l)) (21) 

and therefore 

p>(t)= f drl(r;p(0))p (t-r;\l)). (22) 
Jo 

Differentiation of p = p° + p > gives 

p(t)=p°(t;p(0)) + /(t;p(0))|l)(l|+ f drl(r;p(0))p (t-r;\l)). (23) 



Taking the trace and using trp(t) = 1 gives J(t;p(0)) = A2P22 + A 3 p 33 . From Eq. (|20| ) 
one obtains p°, and inserting this into Eq. (p3|) gives 



p = -j.[Hl cA p - pHU) + (A 2 p 22 + A 3 p 3 3)|l)(l| • (24) 



This is a compact form of the Bloch equations used in Refs. 0. U . 
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In this outline of the quantum jump approach state reductions were used as a tool. 
But it is noteworthy that one can also use the Markov approximation, indicating a close 
connection between the two [26 . 



3. A simple special case: Intermittent probe and ir pulse. 



The quantum jump approach will now be applied to the experimental situation of 
Itano et al. 0. Here one can take A 2 = 0. The simplicity of the mechanism becomes 
particularly clear if the it pulse is switched off while the probe pulse is on. Then Q 2 is 
zero during a probe pulse and Eq. (|16|) reads during this time interval 

/ tt 3 \ 

= -iM . (25) 

V n 3 -iA 3 J 

Note that this annihilates the state |2) and therefore the reduced time development leaves 
1 2) invariant. At the end of a probe pulse one has to wait a short transient time of the order 
of A3 1 for a possible 3-component to decay. This will always be done in the following. 

(i) Effective reduction by a probe pulse 

If the state at the beginning of a probe pulse is 

|V>) = «i|l) + a 2 |2) (26) 

and < r < r p , then evolves, until the emission of the first photon, as 

e-^V/^) = e - MoT |^> = ai e' MoT \l) + a 2 \2) (27) 

since the 7r pulse is assumed to be switched off now. Due to the term A 3 in M the norm 
of the first term of the right hand side decreases exponentially and the first term becomes 
negligible for r large enough. Therefore, if an atom did not emit a photon until the end 
of the probe pulse it will essentially be in the state |2), and this happens with probability 
given by the norm-squared of the r.h.s., i.e. by \ct2\ 2 for large enough r p . 

On the other hand, if an atom does emit one or more photons - this happens with 
probability 1 - | \e~ M ° T 1 | 2 = K| 2 (l - | |e~ M ° T |l) 1 1 2 ) - then right thereafter it is in 
state |1) and will then be pumped between |1) and |3) by the probe pulse, with photon 
emissions. A short time after the end of the probe pulse, |3) decays to |1) due to the 
damping term A 3 . Thus a single atom is projected onto |1) or |2) by the probe pulse with 
probability |«i| 2 and 1 0^2 1 2 if | |e _M ° T |l) 1 1 2 can be neglected. For an ensemble of atoms the 
density matrix becomes diagonal because of the reduction of every single atomic state. 

The preceding analysis is easily made more quantitative as follows. The eigenvalues 
of M are A2 = and 

Ai )3 = \ (A 3 ± ^Al - 4fi 2 ) . (28) 
The first term on the r.h.s. of Eq. (^7]) becomes 

|1> = a x — !-r- {(Mo - A 3 )e- AlT - (M - A^M |1> , (29) 

Ai — A3 1 > 



a ie - Mor ' 
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as immediately checked by explicit differentiation |39[ . For 2f2 3 < A 3 the root in Eq. (|28|) 
is real and Ai and A3 are positive. Therefore in this case the exponential decrease goes at 
least as 



cxp 



T 
1 



A, 



40g 



< exp 



T 



and this becomes exponentially small for a probe pulse of length r p with 

Tp > 2A 3 /tt 2 3 (for 2fi 3 < A 3 ) . 
If AO 2 > Aj the root is imaginary and the decrease goes as 



(30) 



(31) 



exp 



-A 3 
A . 



and Eq. (|31~D is replaced by 



r p > A/A 3 



(for 2Q 3 > A 3 



(32) 



(33) 



This can be combined to 



r p > max I A 3 \ ^3/^3} 



(34) 



For the special case under consideration this is the condition on the length of the probe 
pulse for an effective reduction to |1) and to |2), with probability |ai| 2 and \a 2 \ 2 , respec- 
tively. 



(ii) Population vs. observed photons 

A single atom is projected onto the ground state if it emits several photons during the 
probe pulse. For an ensemble of atoms the number of photons is expected to be a mea- 
sure for the population of level 1. With the quantum jump approach this is easily seen as 
follows. 

The probability for no photon emission until time r, Pq{t; tp), is given by the norm 
squared of the r.h.s. of Eq. (P?[), according to Eq. flT8|), and it approaches |ai| 2 for large 
times. For the subensemble of atoms with emissions the (conditional) probability density 
for the emission of the first photon is therefore 

w 1 (r;^)/\a 1 \ 2 = ~^Po(r; ^)/|«i| 2 (35) 

for ai ^ (for a\ = it is not defined). Since the two terms on the r.h.s. of Eq. 
(f2~7p are orthogonal, the «2 term drops out upon differentiation, and —d/dr Pq(t;iP) is 
proportional to |«i| 2 . Thus for «i 7^ the probability density for the first photon as well 
as the number of photons per atom in this subensemble (i.e. the conditional expectation 
value) is independent of the atomic state at the beginning of the probe pulse. The number, 
N(t;iP), of photons per atom for an unconditioned ensemble with initial state (i.e. 
the usual expectation value) is 

N(t;^) = \ ai \ 2 N(r;l) . (36) 
This expression is now also true for a\ — 0. 
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The result is exact for the case fi 2 = and for all times r within the validity domain 
of the quantum jump approach and Bloch equations. As a consequence, in the case of an 
ensemble realized by a large number of noninteracting atoms without cooperative effects, 
as in Ref . j3| , the number of observed photons per atom is proportional to the population 
of level 1. 



(hi) Effectiveness of state reduction 



Instead of the condition in Eq. ( |34"D for r p one can use iV(r p ; ■?/>), the number of photons 
per atom emitted until time r p , as a more precise measure for the effectiveness of state 
reduction. Eq. (|34]) corresponds to N(t p ; 1) 3> 1, but we will show that also for smaller 
N(r p ;ip) an almost complete state reduction is obtained. 

As pointed out in (i) above, atoms with photon emissions are in |1) if one waits at 
the end of the probe pulse for a short transient time to allow for the decay of level |3). 
But atoms without emissions, however, still contain a part which is not reduced to |2) if 
«i 7^ in Eq. $27\). At the end of the probe pulse and after a short transient time to 
allow for the decay of level |3) this non-reduced component is 

ax (l|e- Mor p|l) |1) . (37) 

The smaller the norm of this, the better the reduction to |2). The norm can be estimated 
by 

\\a x (l|e- A/or p|l) 1 1)| | < ||a ie - MoTp |l)|| . (38) 

For initial state |1) the number of photons per atom until time r p , A(r p ; 1), is in good 
approximation given by the steady state emission rate multiplied by t p , 



N(t p , l) = 4 rjj "a^ ^p- ( 39 ) 

For initial state \ip) = «i|l) +«2|2) we denote the number of photons per atom until time 
r p by N, i.e. N = N(r p ,ip)- Then one can use Eqs. ( |3"6| ) and ( |3"PD to express r p through 
N and lad 
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A 2 4- 90 2 



This can now be inserted into Eq. (RSf) to obtain an estimate of the non-reduced part when 



N photons per atom are emitted. This norm is easily calculated by Eq. (|27j). The norm 
is a function of Q3/A3, N and ai, and one easily shows that for fixed other parameters it 
becomes largest for |ai| = 1. A graphical evaluation gives as upper bound for the norm 
of the non-reduced part 

||ai (l|e- MoTp |l) |1)|| < 1.04 • e - N{T > 1)/2 (41) 

which holds for N > 2 and for all values of f2 3 , A 3 , and a\. For increasing iV the reduction 
thus becomes very effective. 

For particular values of ^3/^3 the bound for the non-reduced part can be substantially 
improved. E.g., for in the vicinity of A%/2 the non- reduced part becomes much smaller 
than the above bound. For very small and very large values of ^3/^3 the reduction is 
somewhat less efficient than for ^3 close to 
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N 


4 


5 


6 


8 


10 


20 


50 


n 3 < a 3 


0.135 


0.082 


0.050 


0.018 


6.7- 10" 3 


4.5 ■ 10" 5 


1.4- 10- 11 


n 3 = a 3 /2 


0.023 


0.006 


0.002 


0.0001 


6.7- 10" 6 


4.0 ■ 10" 12 


2.9 ■ 10" 31 


n 3 = A 3 


0.051 


0.027 


0.015 


0.004 


6.9 • 10" 4 


4.3 ■ 10" 7 


5.1 ■ 10" 17 


= 2A 3 


0.094 


0.065 


0.038 


0.011 


3.4 ■ 10" 3 


1.2 • 10" 5 


5.7- 10" 13 



Table 1: Maximally possible non- reduced part for given number N of observed 
photons per atom and for different values of Q 3 /A 3 . The 7r pulse switched off. 



In Table 1 we have listed the norm of the maximally possible non-reduced part for 
various values of iV and fl 3 /A 3 . The best reduction occurs for Q 3 around A 3 /2, but the 
reduction is also excellent for small and large values of fl 3 if N is larger than 8. 

Summarizing this section, we have shown the following for the case in which the 7r 
pulse is turned off during a probe pulse. 

• The probe pulse provides an effective reduction of the initial state provided its 
duration is much longer than ma.x{l/A 3 , A 3 /Q 3 }, a rather mild restriction. 

• For an (infinite) ensemble the observed number of photons per atom is proportional 
to the population of level 1. 

• Already for small average numbers of emitted photons an almost complete state 
reduction is obtained. 



4. Simultaneous probe and tx pulse 



Now we consider a single probe pulse with the it pulse switched on. At the end of 
a probe pulse, we include a short transient time of the order A 3 X to allow for the decay 
of the auxiliary level 3. Since this short transient time, with the action of the tt pulse, 
is neglected this introduces an error of the order Q 2 /A 3 in the time development of the 
above subensembles. In the following we will assume that this error is much smaller than 
e p . Since e p = A 3 Q 2 /Q 3 this is equivalent to the condition 

e d = tt 2 3 /A 2 3 < 1 . (42) 

In this and the next section we will use the conditions e p , 6d < 1. The condition 6r <C 1 



is then automatically fulfilled f24jj 



The 7r pulse causes a small additional transition between 1 and 2. We will show that, 
as a consequence, an atom with initial state \ifj) = ai|l) + c*2 12) and without photon 
emission until the end of a probe pulse, including the above short transient time, is not 
in the state |2) but in a state |A) which also has a 1-component. On the other hand, 
if an atom emits photons, the last photon may have been emitted some time before the 
end of the probe pulse. Right after the emission the atom is in |1), but until the end of 
the probe pulse a small contribution of state |2) may build up, due to the action of the 
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7r pulse. Thus the atom will not be in |1) as in the ideal projection result. Instead it 
is in a mixed state, denoted by p. Thus, with the 7r pulse switched on, a single probe 
pulse effectively projects onto the state |A) if no photon is emitted and onto p otherwise, 
and this happens with the probability P Q (T p ;^jj) an d 1 — -Po( r p;V0> respectively In the 
following | A) and p will be determined. If e p , <C 1 the differences between |A) and |2), 
p and |1)(1| and Po{t p ; ip) and |ct 2 | 2 are small and the results of Section 3 can be used as 
a good approximation. 

A single atom now evolves with the reduced Hamiltonian H* ed = —ihM of Eq. (|16D , 
with A2 = 0, until the emission of the first photon. The possible pumping between levels 
1 and 2 is reflected by the fact that |2) is no longer annihilated by H^ ed . To calculate 
the time development exp{— i H^ cd t/h} = exp{— Mt} one may proceed in two alternative 
ways if the eigenvalues A; of M, i = 1,2,3, are all distinct. In the first way one determines 
the corresponding eigenvectors |Aj) of M. Since M is non-hermitian these are in general 
nonorthogonal, and therefore one also needs the reciprocal basis {|A*)} with (A*|Aj) = 5ij. 
Then one can write 

e~ Mt = J2 e- Al *|A i )(A i |. (43) 
Alternatively one has, as a generalization of Eq. 



e~ Mt = e~ Xlt ^ — t4tT- — VT + c y clic permutations (44) 
(Ai — A 2 )(Ai — A3) 

which is immediately checked by application to eigenvectors. The case of degenerate 
eigenvalues can be treated by considering appropriate limits of Eq. (pHD . 

Comparing the two equations one sees that |Aj)(A l | coincides with the operator multi- 
plying e~ x,t in Eq. (f44|). Moreover, applying this operator to any vector gives a multiple 
of |Aj), thus automatically yielding the eigenvectors. The eigenvalues are the roots of the 
characteristic equation which, in principle, can be solved in closed form. One easily 
calculates 

/ 1 \ / A 2 (A 2 - \ A 3 ) 
(M-Ai)(M-A 3 ) = |i0 2 (A 2 - \A Z 

V J \ \ 

By the above remarks, this is a multiple of |A 2 ). Similarly 





(M-A0(M-A 3 ) I 1 I =- ^ (M-A!)(M-A 3 ) | | , (45) 

which is also a multiple of |A 2 ). 

For the parameter range of the Introduction good approximations for Aj are 

1 



A 2 = ± A 3 nl/n 2 3 = ±e 2 R A 3 (46) 

where A 2 has been obtained by Newton's method. Note that, for fi 2 <C fi 3 , 

A 2 < Re Ai, 3 . (47) 
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Hence the exponentials exp{— Ai^i} in Eq. ([P]) drop off very rapidly. When calculating 
exp{— Mt p }\i/;) by Eq. (fPJJ one can therefore, as in Section 3, neglect the |Ai) and | A3) 
terms if r p max{A^ , A^/Vt\} = A^/Vtl (because ea <C 1). 



(i) Subensemble without photon emission 

With Eqs. (|44]) to (47) one can now obtain the state of the subensemble of atoms without 
photon emissions until the end HO] of the probe pulse. For initial state 



a 1 |l)+a 2 |2), 



it is described at time r p by 



p 

Mt„l i.\ „-Mt, 



e- MT -{ ai \l) + a 2 \2)} 



e -^r P (aa _ f ) Ji- e l + e 2 | A2 ) ( 4 £ 



where 



-ie p 

r 2 



Moreover, 



Aa>= T. 2 \ I • (49) 

A 2 r p = ^ p tt^«1 (50) 

and hence exp {— A2T P } ~ 1. Similarly for an initial density matrix. The state |A 2 ) in Eq. 
flUf) has a very small third component, and after the probe pulse has been turned off this 
component will, on the time scale of A3 , decay to zero Thus at the end of a probe 
pulse and after this transient decay time the subensemble with no photons is in the - 
normalized - state 

|A> = 2P lj2 |A 2 )/|| • || = -j=L= | 1 I «|2) (51) 

where iPi )2 = |1)(1| + |2)(2| denotes the projector onto the 1-2 subspace and where terms 
of higher orders in e p and e<j have been omitted. The probability for no photon emission 
is 

P (r P ;^) = HJV-^MH 3 

= I ct 2 1 2 + 21m aia^Cp — 7r|a 2 | 2 -7-e p + higher orders in e p , e^. (52) 

For an initial density matrix p m instead of the pure state ip one has to replace |a 2 | 2 by 
p 2 n 2 and «i«2 by p^ 2 . 

Thus, to good approximation, the probability for no photon emission is proportional 
to the population of level 2 and the atoms with no emissions are approximately in the 
state 1 2). For the parameters [44] of the experiment || the corrections in Eq. (|52|) are 
less than 4 • 10" 4 . 
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(ii) Subensemble with photon emissions 



We will now calculate the density matrix p for the subensemble with photon emissions. 
One can employ a systematic expansion in powers of Sl 2? including second order. However, 
the following more physically motivated procedure is simpler and yields the same results. 
At the beginning of the probe pulse we assume the ensemble to be in the pure state 
for a density matrix the treatment is similar. The complete ensemble at time r p after 
the beginning of the probe pulse can be thought of as consisting of two subensembles of 
atoms with and without photon emissions. The latter is described by Eq. (|20|) , 

p (r p ,^^e~ M ^)(i;\e- M ^ , (53) 

with its relative weight given by Eq. (|18|), 

trp°(r p ;^) = P (r p , ^) . (54) 

According to Eq. (|22"D the former subensemble is described by 

p>(r p ; V) = £ dr I(r; ^)p°(r p - r; 1) . (55) 

The complete density matrix is p > + p°, and therefore 

trp> = 1 - P . (56) 

We now determine p> 2 and p% 2 at time r p . For initial state |1) the no-photon probability 
decreases rapidly, and therefore p°(r p — r; 1) contributes essentially only in the vicinity of 
r « r p (cf. Eqs. fl2"9|) to (|3T|)). Because of this one can replace I(t;i[)) by I(t p ;i/;) which 
is approximately equal to the fraction 1 — Pq(t p ; ip) °f emitting atoms times I(t p ; 1). The 
latter practically equals the stationary emission rate for the three-level system which is, 
up to terms of order in the stationary rate from Eq. (|39D for the two-level system. For 
the calculation of the 12 and 22 component these corrections in fl| can be omitted since 
p (r p — r; l)i 2 and p°(r p — r; 1) 2 2 are themselves of order Q 2 and f2| and much smaller 
than 1. Thus we obtain 

P > (r P ;^) 1 2/22 = (1 - P ^ r p^)) A 2 3 + 2Ql £ rfrp ° (r; 1)12/22 • (57) 

After the end of the probe pulse, any population of level 3 will rapidly decay to level 1 
in a transient time of order A% . We denote the resulting - normalized - density matrix 
of the subensemble of atoms with emissions by p [fi~5|| . The normalization factor of p > is 



1 — Pq, and according to Eq. ( |56"D we obtain from Eq. (|5 

A 3 n 2 3 



P(r p ; V»)i2/22 = , 2 1 002 / " dr P°( T ; 1)12/22 • (5? 



For this we note that, by Eq. ( p3|) , 

p? 2 (r;l) = (l|e-^|l>(l|e-^|2> 

p° 2 (r;l) = |(2|e-^|l)| 2 - (59) 
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Using Eq. (f44|) for e Mr an elementary calculation yields 

P22 = tt^t- e P + higher orders in e p , e d (60) 

P12 = + higher orders in e p , e d (61) 

pn = 1 - P22 , P21 = - P12 (62) 

and its 13, 23 and 33 components vanish, 

P13 = P31 = P23 = P32 = P33 = 0. (63) 

We conclude that after a short transient time at the end of the probe pulse the 
subensemble with photon emissions is described by the above normalized state p, which 
is independent of the initial state and the relative size of the subensemble is given by 

l-P (r P ;V)- 

We note that p 22 arid p\2 are indeed very small. For the parameters of the experiment 



H one has p 22 < 1.2 • 10~ 5 and \p 12 \ < 4.1 • 10" 



(iii) Level population after a probe pulse 

We denote by p^ the density matrix of the complete atomic ensemble after a short 
transient time at the end of a probe pulse. By the preceding results it is given by 

p( p ) = P (Tp;p in )|A>(A| + (1 - P (r P ;p in ))p (64) 

where p m is the density matrix at the beginning of the probe pulse. From Eqs. (0), (|51~1) 
and (|60D one immediately obtains for the population of level 2 after the probe pulse 

= p 2 \ + 2e p Imp* + vr^e p (l - 2p 2 n 2 ) . (65) 

The first term is the projection-postulate result for an ideal measurement. For n probe 
pulses, P12 is of the order sin | ( ^ — ^-J , and if n is as in the experiment || the second 
term is larger than, or comparable to, the last term. In the corresponding Eq. (16) of 
Ref. H this important term is missing, due the approximation used there f22|, and the 
correction appearing there is equivalent to 0(e p 27rr p /T 7r ) in our notation. 



Summarizing this section, we have shown that, for the parameters e p and much less 
than 1, a probe pulse acts as an effective state reduction also in the presence of the 7r 
pulse. The reductions are to p and |A), corresponding to subensembles with and without 
emissions, respectively. One has p ~ |1)(1| and |A) ~ |2). The corrections have been 
explicitly calculated in terms of the above parameters [f5|. 



5. Level population after n probe pulses 



The preceding results show that, after a probe pulse and a short transient time, the 
ensemble consists of two subensembles, one in the state p and the other in |A), corre- 
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sponding to atoms with and without photon emissions. The small difference from the 
projection result can have a cumulative effect for the density matrix after n probe pulses. 
After a probe pulse, the density matrix for the complete ensemble is of the form 

p = ap + (3 |A)(A] (66) 

with a + /3 = 1. Until the beginning of the next probe pulse at a time T n /n — r p later the 
time development is given by the ir pulse only, i.e. in matrix notation and in the |1) — |2) 
subspace by 

UAt,0)=( C ° S ^ f --ini^t\ (67) 
v ' \ -«sm | \l 2 t cos | ll 2 t J y ' 



For t = — — r p we define U n (t, 0) = £/ n . We will now determine the density matrix after 
the A;-th probe pulse. To this end we put 

p = tr{P h2 e- Mr *U n pUle- MtT vlP lt2 } 

q = \\P lt2 e~ M ^U n \\)\\\ (68) 

where iPi j2 is the projector onto the 1 — 2 subspace. Physically, p is the probability of 
finding no photons after the next probe pulse [||| if one had started with p at the end of 
the preceding probe pulse. Similarly for q and |A). 
With the abbreviation 

s n = sin 7T — 

\n In 

c n = cos Ti ( -7-) (69) 



one finds by a straightforward calculation from Eq. ( |6"TD for U n and Eq. for e~ MTp 

1 1 r 

p = -(I - c n ) + s n e p + p 22 c n - ipi 2 s n - -tt^-(1 - c n ) e p 

= t;(1 - c n ) + 2s„,e p + 7r^-c„e p - t^|t(1 - c„) e p 

1 1 r 

9 = 2( 1 + c ") - 2s «e P - 2 7r 2^"( 1 + c n) e P ( 70 ) 

where higher orders in e p and ea have been omitted. 

Now suppose that after the (k — l)-st probe pulse the density matrix p is given by Eq. 
(H|) with a = a(k — 1) and [3 = (3{k — 1). Then after the k-th probe pulse the relative 
size of the no-photon subensemble is given by 

p(h) = p a(k - 1) + q /3(k - 1) 

= p+(q-p) (3(k-l) (71) 
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where a = 1 — (3 has been used. The solution of this recurrence relation is 



a{k) 



P 



(q-p) 



k-l 



I- {q-p) 
l-P(k) . 



+ {q~pf- 1 m 



(72) 



According to Eq. (|66|) is the probability to find no photon during the first probe 
pulse [0 . For the initial condition that all atoms are prepared in the ground state at the 
beginning of the experiment, as in the experiment of Ref. ||, f3{l) is given by 



P h2 e- MT *U n \l) 



-(1 - c n ) + s n e p - ^^(1 - c„) e p . 



At the end of the ir pulse, i.e. after the n-th probe pulse the density matrix is 

p(T„) = a(n) p + (3(n) |A)(A| . 
The populations of levels 2 is then 

p 22 (T n ) = (2\p(T v )\2) = a(n)fa + /3(n) |(2|A)| 2 . 



(73) 



(74) 



(75) 



Using Eqs. 
obtains 



and Eqs. (|70|)-(|75D and omitting higher orders in e p and e d one 
1 



;i - cl) + (2n - l)s n c n -\ + Trn-f c n n e 



T 



1 - / o 22 (T 7r ) . 



(76) 



Omitting also the e p terms one obtains for the populations of levels 1 and 2 the approxi- 
mate results 



P22(T 7r 



Pu(T n ) 



2 (1_C " } 



= + <) 



1 - COS n 71 



n 71 



1 + COS" 7T ( - - ^ 



(77) 



For n ideal measurements with the projection postulate the result would be |3], 



P22(T 7T ) 



Pn(T n ) 



1 

2 



1 - cos r ' 
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1 + cos" 



7T 



(78 



This differs from the result in Eq. ( |77| ) only by the term r p /T n in the cosine, and therefore 
Eq. (|77|) can be obtained from the ideal projection result by neglecting the action of the 
7r pulse during a probe pulse, thus replacing the time AT = T^/n between measurements 
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36056 
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335 


8 


0.23460 


0.20857 


0.20998 


0.20993 


0.194 


16 


0.13343 


0.10029 


0.10215 


0.10212 


0.103 


32 


0.07156 


0.03642 


0.03841 


0.03840 


0.013 


64 


0.00371 


0.00613 


0.00789 


0.00789 


-0.006 



Table 2: Predicted and observed population of level 2 at the end of the it pulse for 
n probe pulses of length r p . 



by the effective time AT = T n /n — r p . This can be also understood directly quite easily 
as follows. Atoms which emit photons during a probe pulse flip repeatedly between levels 
1 and 3, and so the 7r pulse acts less effectively them. Moreover, right after emission of a 
photon an atom is in the ground state and since the action of the ir pulse is of cosine form 
it is small for small times. Similarly, as shown in Section 4, the atoms without emissions 
rapidly approach the state |A) = |2) so that one again has a small action for small times. 
Therefore the action of the it pulse is greatly inhibited during a probe pulse . 

It is evident from Eq. ([77]) that this approximation yields the same result as if one 
had switched off the ti pulse during the probe pulse and then uses the results of Section 
3. 

The corrections of Eq. ( [76] ) to the approximate values for pu(T n ) in Eq. ([H]) are 
small for e p « 1. Moreover, it is straightforward to show that in the parameter range 
considered here the correction to p22{T n ) is positive and increases with n, as long as r p 
is not too close to T T /n. This is borne out in Table 2 where predicted and observed 
values of p22(Tr) for the parameters of the experiment of Ref. || are shown. The second 
column is based on the projection postulate for n ideal measurements. The third column 
is based on Eq. ([77]) or, alternatively, on n ideal measurements with ensuing switch-off of 
the n pulse for r p seconds. The agreement between the quantum jump result in Eq. (|76| ) 
and the numerical solution of the three-level Bloch equations of Eq. fl24|) in column 5 is 
apparent. The projection postulate with AT modified to AT = T w /n — r p also gives very 
good results. The experiment deals with a system where additional energy levels may 
make minor contributions and this may explain the deviations from the experimental 
results in the last column of Table 2. 



6. Conclusions. Does the Zeno effect exist? 



We have investigated the so-called quantum Zeno effect for an ensemble of atomic 
three-level systems as that of the experiment in Ref. J| . There has been some controversy 
about the interpretation of that experiment as to whether it provides an experimental 
proof of that effect. The Zeno effect is a theoretical prediction for the behavior of a system 
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under rapidly repeated measurements with ensuing state reductions according to the 
projection postulate. As explained in the Introduction, the controversy has mainly focused 
on two points, namely (i) whether some of the "measurements" in Ref. [|3] should not 
rather be included as "interactions" in the Hamiltonian, and (ii) whether the projection 
postulate is appropriate at all. 

Now, a measurement on a microscopic system is indeed a complicated interaction with 
another system, ultimately at the macroscopic level. Not all conceivable measurements 
conform to the idealized case considered by von Neumann and Liiders || where each 
measurement is associated with an operator, A say, such that the mean value, (A), is 
given by the expectation value and the mean-square deviation by the expectation of 
(A — (A)) 2 . As a simple example one may consider a system of photons in a cavity. The 
outcome of the measurement consists of two possible pointer readings, and 1 say. The 
actual measurement is performed by a "black box" which has been so constructed that 
an auxiliary two-level atom in some initial state is passed through the cavity and then 
it is determined whether the atom is in the ground or excited state, yielding the pointer 
readings or 1, respectively. Considered as a measurement performed by the black box 
on the photon system, it is obviously not of the above kind. There is indeed an operator 
A for the photon system, whose expectation value gives the mean, but it has eigenvalues 
different from and 1, and right after a measurement the photon is not in an eigenstate of 
A. As a consequence, the mean-square deviation is not of the above form. By combining 
part of the black box, namely the atom, with the photons to a larger system one may 
possibly retrieve or come close to the situation considered by von Neumann, depending 
on how the measurement on the atom is actually performed. 

This example shows that one may have to study the measurement at hand more 
closely. This is what we have done with the atomic Zeno experiment by means of the 
quantum jump approach ]26 , [27j, which is essentially equivalent to the Monte-Carlo wave 
function approach and to the quantum trajectory approach [P^| . The advantage of 
this approach is that it allows a physically intuitive as well as analytic treatment of the 
problem. Part of the measurement interaction - the laser pulses - has been incorporated in 
the Hamiltonian; this corresponds in the above example to combining the photon system 
and atom to a larger system and is related to shifting the Heisenberg "cut" |^ . 

Our analysis has shown that, to a fair accuracy and within the parameter regimes 
considered in this paper, a short probe pulse can indeed be viewed as performing a mea- 
surement of levels 1 or 2 of the atom, with ensuing state reduction, as given by the 
projection postulate for an ideal measurement. However, since these "measuring" pulses 
have been modeled very accurately in the Hamiltonian we were able to show that the 
more realistic case is also slightly more complicated, giving rise to corrections to the 
idealized case. These corrections were explicitly calculated, and they have a cumulative 
effect on the density matrix when the number of probe pulses is increased - i.e. just for 
the interesting case. 

What then remains of the Zeno effect? Does it exist at all? In our opinion the 
answer depends on one's point of view. If one takes the view that for example the probe 
pulses ("measuring pulses") have nothing to do with a measurement but just lead to 
additional terms in the Hamiltonian, then any change in the temporal development is 
not surprising and may simply be attributed to these additional interaction terms. The 
other - possibly more fruitful - point of view is that these pulses approximately realize 
measurements with state reductions, and then one immediately has simple predictions for 
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the approximate behavior of the system and arrives at the impediment and slow-down 
of the time evolution without complicated calculation. Finer details require of course a 
finer analysis, as performed in the previous sections. An actual freezing of the state does 
not seem possible since all realistic measurements take a finite time. In the present case 
this hinges on the corrections and on the finite duration of the probe pulse (including the 
transient decay time of order A3 1 ). 

Our analysis may possibly shed some light on the use of the projection postulate in 
quantum optics in general, not only in connection with the Zeno effect. It seems that quite 
often the projection postulate is a useful tool which can give quick and fairly accurate 
answers. The accuracy depends on how far the particular realistic measurement differs 
from an ideal measurement as considered in orthodox quantum mechanics, and corrections 
may have to be taken into account. The idealization of realistic measurements and the 
projection postulate may often be very useful. Over-idealization, however, should be 
avoided since any idea, when carried to extremes, easily reduces itself to absurdity. 

Acknowledgments. One of us (G.C.H) is indebted to the late Gerhart Liiders (1920 
- 1995), who also discovered the TCP Theorem, for stimulating discussions on his formu- 
lation of the projection postulate. 
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After completion of this paper we received a preprint on the QZE by W.L. Power 
and P.L. Knight with numerical simulations of quantum trajectories. 



Another effect intimately related to the Zeno effect is the so-called inverse quantum 
Zeno effect proposed by V. Aharonov and M. Vardi, Phys. Rev. D 21, 2235 (1980) 
according to which a system can be made to follow an arbitrarily chosen trajectory by 
performing a suitable sequence of measurements. T.P. Altenmuller and A. Schenzle, 
Phys. Rev. A 48, 70 (1993) have discussed this inverse effect for a 2 level system. 

Ref. ]TJ| also uses extensive numerical simulations of quantum trajectories to gain 



insight into the evolution of the atomic-level population when the system is contin- 
uously observed. We, on the other hand, use quantum trajectories in an analytic 
way. 

From this one can easily calculate the populations of levels 1 and 2 after a single 
probe pulse, and they turn out to differ from Eqs. (11) - (16) of Ref. || by an 
additional correction to the projection-postulate result (see our Eq. fl65]) below). This 
is due to the approximation used in Ref. f| for the solution of the Bloch equations. 
This leads to a rate equation for p 22 in their Eq. (9) in which p 12 does not appear, 
and it introduces an error of the order of the missing term as pointed out by D.G. 
Sondermann (private communication). 

This is seen as follows. The second condition means that on average several photons 
should be produced during a probe pulse. Estimating the number of photons by the 



steady-state solution of the Bloch equations [32] for the 1 «-> 3 transition, driven by 



the probe pulse and with no population in level 2, gives the condition 



where Sl 3 is the Rabi frequency of the probe pulse. This last condition now implies 



A 3 2 



3 



and hence, by Eq. (|IJ), 



7T A3 2 

T = — S> — - + — 

n n 2 y> n§ + a, ■ 



This implies the first part of Eq. (|3|). The last equation can also be written as 

n 3 \n 3 a 3 J 

Using the inequality x + 2/x > 2v^2 for x > this gives the condition on the Rabi 
frequencies in Eq. (|3]). This condition is also directly obtained by noting that level 
1 has to be pumped at least once to level 3, and hence r p > n/Q 3 . Eq. (|l|) then 
implies cr «C 1. This condition also automatically ensures the obvious requirement 
that during a probe pulse the ir pulse should have negligible effect. 
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Figure 1: V system with (meta-) stable level 2 and Einstein coefficient A 3 for level 3. Q2 
and fi 3 are the Rabi frequencies of the two lasers. 
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Figure 2: Probe pulses and ix pulse 
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Figure 3: 7r pulse switched off during probe pulse 
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